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In present report we review some preliminary results of investigation of higher orders of mean-field expansion 
for Nambu— Jona-Lasinio model. We discuss first results of investigation of next-to-next-to-leading order of 
mean-field expansion equations for four-quark and three-quark Green functions. We have considered equations 
for Green functions of Nambu-Jona-Lasinio model in mean-field expansion up to third order. 
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1. INTRODUCTION AND SUMMARY 



Multi-particle equations are a traditional basis for 
a quantum-field description of bound states in parti- 
cle physics. A well-known example of these equations 
is the two-particle Bethe-Salpeter equation (BSE) for 
the two-particle amplitude and for the two-particle 
bound state Q.The multi-particle (three or more par- 
ticle) generalizations of the BSE have been also stud- 
ied. A straightforward generalization of two-particle 
BSE has been intensively studied in sixties-seventies 
of last century. A best exposition of these studies can 
be found in the work of Huang and Weldon ^ . 

An essential imperfection of the original Bcthc- 
Salpeter approach to multi-particle equations was a 
full disconnection of the approach with the field- 
theoretical equations for Green functions (which are 
known as Schwinger-Dayson equations (SDE)). This 
imperfection has been eliminated by Dahmen and 
Jona-Lasinio, which had included the BSE to the field- 
theoretical Lagrangian formalism with the consider- 
ation of functional Legendre transformation with re- 
spect to bilocal source of fields . Then this approach 
has been generalized for multi-particle equations with 
consideration of Legendre transformation with respect 
to multi-local sources Q. 

However, these theoretical constructions had not 
solved the principal dynamical problem of quantita- 
tive description of real bound states (nucleons, mesons 
etc.). A solution of BSE- type equations has been 
founded as a very complicated mathematical tool even 
for simple dynamical model. There is a main rea- 
son of a comparatively small popularity of the method 
of multi-particle BSE-type equations among the theo- 
rists. Much more popular approach to the problem of 
hadronization in QCD is based on the 't Hooft's con- 
jecture that QCD can be regarded as an effective the- 
ory of mesons and glueballs Subsequently, it was 



shown by Witten that the baryons could be viewed as 
the solitons of the meson theory 01 . Futher develop- 
ment of these ideas has been successful and has leaded 
to theprediction of pentaquark states in baryon spec- 
trum [3. 

Nevertheless, the investigations of multi-quark 
equations are of significant interest due to the much 
less model assumptions in this approach in comparison 
with the chiral-soliton models. The solutions of multi- 
quark equations will provide us almost exhaustive in- 
formation about the structure of hadrons. There is 
basic motivation of present work. 

We shall investigate Nambu-Jona-Lasinio (NJL) 
model with quark content which is one of the most 
successful effective models of QCD in the non- 
perturbative region (for review see over- 
whelming majority of the investigations, the NJL 
model has been considered in the mean-field approx- 
imation or in the leading order of l/uc- expansion. 
However, a number of perspective physical applica- 
tions of NJL model is connected with multi-quark 
functions (for example: meson decays, pion-pion scat- 
tering, baryons, pentaquarks etc.). These multi-quark 
functions arise in higher orders of mean-field expan- 
sion(MFE) for NJL model. To formulate MFE we 
have used an iteration scheme of solution of SDE with 
fermion bilocal source ^3 ■ 

We have considered equations for Green functions of 
NJL model in MFE up to third order. The leading ap- 
proximation and next-to-leading order (NLO) of MFE 
maintain equations for the quark propagator and the 
two-quark function and also the NLO correction to the 
quark propagator. The next-to-next-to-leading order 
(NNLO) of MFE maintains the equations for four- 
quark and three-quark functions, and next-to-next- 
to-next-to-leading order(NNNLO) of MFE maintains 
the equations for six-quark and five-quark functions. 
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2. MEAN-FIELD EXPANSION IN 
BILOCAL-SOURCE FORMALISM. 
LEADING ORDER AND NLO 

We consider NJL model with the Lagrangian 



The Lagrangian is invariant under transformations 
of chiral group SUv{2) x SUa{'2), which correspond to 
u-d quark sector. A generating functional of Green 
functions (vacuum expectation values of T-products 
of fields) can be represented as the functional integral 
with bilocal source: 



Girj) = J D(V',i/')expi| J dxC 

dxdyip{y)Tj{y,x)ip{x)Y 

Here ri{y, x) is the bilocal source of the quark field. 

The n-th functional derivative of G over source 77 is 
the n-particle (2n-point) Green function: 



ri=0 



< OT 
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Translational invariance of the functional-integration 
measure gives us the functional-differential SDE for 
the generating functional G : 
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We shall solve this equation employing the method 
which proposed in work . A leading approximation 
is the functional 



= exp 



The leading approximation generates the linear itera- 
tion scheme: 

G = G(°) + G(i) + ■ • • + G(") + • • ■ , 

Functional G*^"^ is 

Q{n) „ p(")(;;(o)^ 

where P*^"^ is a polynomial of 2n -th order over the 
bilocal source 77 . 

The unique connected Green function of the lead- 
ing approximation is the quark propagator. Other 
connected Green functions appear in the following it- 
eration steps. The quark propagator in the chiral limit 
is 

where m is the dynamical quark mass, which is a so- 
lution of gap equation. 
A first-order functional is 

G« = |iTr(5W*r;2) + Tr(5« * 77) }g(°). 

The iteration-scheme equations give us the equation 
for two-particle function ^2^^: 

-Sq{x - y')S(){x' - y) 



(1) / 2: y \ _ 



x' y' 



ig dxi^{Sa{x-xi)Sa{xi~y))tr 



5f) 


f Xi Xi\ 




U' y' ). 



' ^ ^ X y 

(1) 



-iSo{x-xi)j5T''T'''Soixi-y))tr 

and the NLO correction to quark propagator S'^^ 

s<"(=^-rt = is/<i«s™(---.){s5"(^; I 

+ig [ dxiS^°H X — Xi )S^°\xi-y)trS'-^\0). 



The graphical representations of these equations see 
on Figs. 1 and 2., where the graphical notations of 
Fig. 3 are used. 

These equations can be easily solved, and the so- 
lutions contain singlet scalar quark-antiquark bound 
state with mass 2m (sigma-meson) and massless (in 
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Figure 1: The equation for NLO two-quark function. 



Figure 4: The equation for NNLO four-quark function. 



Figure 2: The equation for NLO correction to quark 
propagator. 




the chiral limit) pseudoscalar bound states (pions). 
To describe the solution of the NLO equation for two- 
quark function and for future purposes we introduce 
the composite meson propagators by following way: 
a) Let us define scalar-scalar function 



Sa{x — x') = tr 



{I) I X X 



ipip{x)ip'ip{x') > 

From the equation (1) for two-quark function we ob- 
tain (in momentum space) 



Figure 5: The equation for NNLO three-quark function. 



- T - T 

~< V'75Y'0(a;)'075YV'(a;') > (5) 

From the equation (1) for two-quark function we ob- 
tain (in momentum space): 



5f(p) = -l(<5''^-zAf(p)). 
Here we define the pion propagator 



(6) 



^9 



(3) 



Here we define the following function, which we call 
cr— meson propagator 



Af(p) 

where Z^ip) = AW . 



S''''Z{p) 



(7) 



where 



A.(p) 



zAp) 



Z{P) 



and 



hip) ^ dq 



4m^ — p^ ' 
/o(4r77?) 

1 



(4) 



(m^ — {p + qY){m'^ — q^) 
b) Pseudoscalar-pseudoscalar function is defined as 



S-f (.T -x')=tr 



S. 



2 \ x' X' l^'Y^^Y 
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Figure 3: Diagram rules. 



3. NNLO OF MEAN-FIELD EXPANSION 

NNLO(second-order)generating functional is 



1 
3! 



The equations for four-quark and three-quark func- 



tions see on Figs. 4 and 5. 

(2) 

The equations for the four-quark function S\ and 

(2) 

for the three-quark functions 6*3 are new, and the 

(2) 

equations for two-particle function 5*2 and propa- 
gator 5'^2) have the same form as the corresponding 
NLO equation except of the inhomogeneous terms. 
For NNLO equations these terms contain four-quark 
function s'^^ and three-quark S'g^'' function. 

The equation for the four-quark function has a sim- 
ple exact solution which is the product of NLO two- 
quark functions (see Fig. 6). As it seen from this solu- 
tion, the TTTT— scattering in NJL model is suppressed, 
i.e. in the NNLO of MFE this scattering is absent, 
and it can be arise in the NNNLO (third order) only. 
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Figure 6: The solution of equation for four-quark 
function. 
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Figure 7: The connected part of sigma-pion-pion- vertex. 



4. VERTEX amv 

The existence of above exact solution for the four- 
quark function gives us a possibility to obtain a closed 
equation for the three-quark function. As a first step 
in an investigation of this equation we shall solve a 
problem of definition of crTTTr— vertex with composite 
sigma-meson and pions. Let us introduce a function 



W^U^x'x") EE tr 



.(2) 



75 ^75 IT 
x" x" ' ^ ^ 



- T - T 



and define: 

a) scalar vertex 



V„(xx'x")=tr 



So{x - x')S^^'' \ 5/ 5' 



2inc / dxiVs{xx'xi)A„{xi — x"). (8) 



Here: vs{xx'x") = tr4So{x-x')SQ{x' -x")Sq{x" -x)] 
is the triangle diagram, 
b) pseudoscalar vertex 



V^'ixx'x") = tr 



So{x-x')-f5 — S. 



2 ~2 I x" x" j ^5 2 



= 2mc j dxivp{xx'xi)/^l\xi - x"). (9) 

Here: vp{xx'x") — tra[Sa{x — x')^zSq(x' — 
x")-i,S^{x" - x))]. 

With definitions (2)-(9) we obtain for vertex func- 
tion W"^'^ the following equation: 

W:l^{xx'x") = Wtixx'x") 



+2igncJ dxilsix - xi)W^^^{xix'x"), 



where ls{x) = tra[So{x)So{~x)] is the scalar quark 
loop. Inhomogeneous term Wq'^ is: 

WS^xx'x") = Vfixx'x") + Vf{xx"x') 



+Aig I dxiV^'''{xxix')S^'\xi ~ x") 



+Aig I dxiV^'''ixxix")S^'''{xi - x') 



ig / dxi{V^{xxiXi) - AV^"''{xxiXi))Sf{x' - x") 



Using definitions (2)- (9) we have: 

[w;}^ixx'x")Y°'' 



= -2nc J dxidx2Vp{xxiX2)[A'^'''{x2-x')A'^''''{xi-x") 



+Al''"'{x2 - x")A-^'''{xi - x% 

The equation for W""^ can be easy solved in the 
momentum space and and a solution is 

W^'^^ippY) = tA,ip)W^\pp'p") 

where p is cr— mesons momenta, and p' , p" are pion 
momenta: p = p' + p" . The connected part of W"^^ 
is an decay amplitude a — > tttt. It is has a following 
form: 



On 

= -^A,{p)[vp{pp'p")+vp{pp"p')]AT' [p')Al-\p"), 

(10) 

(See also Fig.7.). 



5. NNNLO OF MEAN-FIELD EXPANSION 

The third-order generating functional is 
G(^)W = {^Tr(4^U,«)+iTr(^,f)*,^) 



IPM-LHP06-19May 



IPM School & Conference on Lepton and Hadron Physics, Tehran, May 15-20 2006 



5 




Figure 8: The equation for six-quark function. 




Figure 9: The equation for five-quark function. 
+ iTr(5f *,^)++iTr(5f 

The graphical representations of equations for six- 
quark functions and for five-quark function see on 
Figs. 8 and 9. The equations for the six-quark 
function and for the five-quark function in our iter- 
ation scheme are new, and the equations for four- 

(3) . (3) 

quark function S\ , three-quark function S\ , two- 
quark function 5*2^' and quark propagator S^^^ have 
the same form as the second-order equations except of 
the inhomogeneous term, which contains the six-quark 
function and the five-quark function. 
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